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The s = l/2 chain-boundary excitations occurring in the Haldane phase of s = 1 antiferromagnetic 
spin chains are investigated. The bilinear-biquadratic hamiltonian is used to study these excitations 
as a function of the strength of the biquadratic term, /3, between — 1 < /3 < 1. At the AKLT point, 
ft = — 1/3 we show explicitly that these excitations are localized at the boundaries of the chain on 
a length scale equal to the correlation length £ = l/ln3, and that the on-site magnetization for the 
first site is < Sf >= 2/3. Applying the density matrix renormalization group we show that the 
chain-boundary excitations remain localized at the boundaries for — 1 < /3 < 1. As the two critical 
points P = ±1 are approached the size of the s = 1/2 objects diverges and their amplitude vanishes. 



Many different compounds exist which are believed to 
be quasi one-dimensional s = 1 antiferromagnetic sys- 
tems. Among these are CsNiCL;, which yielded the first 
experimental confirmation 1 of the Haldane gap 2 , as well 
as Ni(C 2 H 8 N 2 ) 2 N0 2 (C104) (NENP) 3 ' 4 which has been 
studied extensively and Y 2 BaNiOs 5 , a compound sim- 
ilar to the perovskites, which has garnered considerable 
attention recently due to the relative ease with which it 
is carrier doped. These compounds are thought to be in 
the so-called Haldane phase at low-temperatures. This 
phase is characterized by a disordered singlet ground- 
state with a gap to the lowest lying excitations in the 
thermodynamic limit. In general, antiferromagnetically 
coupled linear s = 1 systems are believed to be in the 
Haldane phase if the hamiltonian describing the interac- 
tion between the spins is sufficiently close to the stan- 
dard nearest neighbor Hciscnbcrg model. In particular it 
is known that if a biquadratic term linear in (3 is added: 

L-l 

H = J Y, [S,:-S 4+1 -/3(S,-S l+1 ) 2 ] , (1) 

i=i 

the system remains in the Haldane phase for — 1 < (3 < 1. 
At the Affleck-Kennedy-Lieb-Tasaki 8 (AKLT) point, (3 = 
— |, this model is partly solvable and the exact ground- 
state corresponds to a valence bond solid (VBS) of adja- 
cent dimers 8 (Fig. 1) and throughout the Haldane phase 
the ground-state is presumably closely related to the VBS 
state. If each s — 1 spin is viewed as two s = 1/2 spin 
in the symmetric triplet state then the VBS state cor- 
responds to combining each of the two s = 1/2 into a 
singlet with spins on the adjacent sites. 

As was noted by Kennedy 9 , it is only when periodic 
boundary conditions are imposed that the ground-state 
can be considered as a non-degenerate singlet. If the 
chain is broken (i.e. open boundary conditions) then the 
ground-state is actually four-fold degenerate in the ther- 
modynamic limit. One way to understand this four- fold 
degeneracy is to think of the open chain as a periodic 
chain with one coupling removed. At the AKLT point 
this corresponds to breaking one of the dimers, result- 
ing in two free s = 1/2 objects and a ground-state that 




figure 1: Schematic representation of the VBS state for 
open boundary conditions. Large circles correspond to 
spin-1 sites and the smaller one to spin-1/2 states. The 
s = 1/2 excitations at the boundaries are indicated by 
the bold arrows. 

is four-fold degenerate in the thermodynamic limit, as 
shown in Fig. 1. Hence, it was suggested 10 that such 
s = 1/2 are real and experimentally observable. For a 
finite open chain the s = 1/2 excitations couple and an 
exponentially low- lying triplet, above the singlet ground- 
state, is found in the Haldane gap 9 . In the Haldane phase 
these s = 1/2 objects can be shown to be located at the 
boundaries of the open chain and we shall refer to them 
as chain-boundary excitations. 

The chain-boundary excitations may seem like a cu- 
riosity, but it is important to realize that the Haldane 
systems are gapped and hence, in the presence of impu- 
rities, the low-temperature physics will be determined by 
the interaction between these s = 1/2 objects. In partic- 
ular, there is considerable current interest in disordered 
s = 1 chains 11-14 which in many cases are treated in 
terms of interacting s = 1/2 objects. Hence, considerable 
effort has been devoted to the experimental verification of 
the above considerations. Electron-spin-resonance (ESR) 
studies on Cu doped NENP 10,15 have confirmed the exis- 
tence of these s = 1/2 chain-boundary excitations. Sub- 
sequent studies of NENP doped with non-magnetic im- 
purities Zn, Cd and Hg also showed evidence for s = 1/2 
excitations 16 . More recent specific heat measurements 17 
on doped Y 2 BaNiOs has been interpreted as low- lying 
s — 1 excitations, but a more detailed analysis showed 
that an interpretation in terms of s = 1/2 excitations 
also is possible 18 . 

In light of the above remarks it seems important to 
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understand how the chain-boundary excitations occur 
throughout the Haldane phase, a point we systemati- 
cally investigate in the present paper. In order to include 
exact results we have studied the bilinear-biquadratic 
model, Eq. (1), with the physical systems correspond- 
ing to (3 — 0, and we take this model as representa- 
tive of the Haldane phase. At the AKLT point the 
s = 1/2 excitations can be studied using the exact re- 
sults for the VBS state as we detail below, facilitating 
the analysis. Exponentially localized s — 1/2 objects are 
shown to occur at each end of an open chain. Using the 
density matrix renormalization group (DMRG) we show 
that these s = 1/2 chain-boundary excitations remain 
well-defined in the Haldane phase as (3 is varied between 
— 1 < /? < 1, disappearing as the critical points (3 = ±1 
are approached. 

We now briefly discuss the physics of the Hamilto- 
nian Eq. (I). Two critical points exist. The Uimin- 
Lai-Sutherland (ULS) point 19 , described by a SU(3)fc = i 
Wess-Zumino-Witten (WZW) model, at (3 = —I, where 
the system displays a Berezinskii-Kosterlitz-Thouless 
phase transition with an exponentially diverging corre- 
lation length into a massless phase 20 for (3 < — 1, and 
the Takhtajan-Babujian (TB) point 21 , described by a 
SU(2)fc =2 WZW model, at (3 = 1 where a second order 
phase transition to a dimerized phase occurs. At the ULS 
point the system has gaplcss modes at k = 0, ±2tt/3 in 
contrast to the Heisenberg point, (3 = 0, where the Hal- 
dane gap occurs at k = it. Consequently, as the ULS 
point is approached, a disorder point, which has been 
identified with the AKLT point 22 , occurs where the sys- 
tem develops short-range incommensurate correlations. 
The associated Lifshitz point, where the peak in the 
structure factor moves away from k = tt, has been nu- 
merically estimated to occur at (3 ~ — 0.438 22,23 . At the 
TB point gapless modes occur at k = 0, ir and hence, as 
(3 approaches 1 the peak in the structure factor remains 
at 7r. 

Exact Results: The VBS ground-state can be writ- 
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Here |+ >i,\0 >i,\— >i corresponds to the three states 
of the s = I on site i. If periodic boundary conditions 
are considered the trace should be taken. Using the 
above wave-function it is now straightforward to eval- 
uate < Sf >vbs, the magnetization on the first site of the 
open chain, as a function of L. One finds for the state 
with total magnetization Sj. = I: 



< Sf > v bs (L) = 



2-3" 



I - 3- 



L even 



(3) 



for even L, including all finite size corrections. As is 
evident from this expression this amplitude rapidly ap- 
proaches 2/3 as L — > oo. For odd L one finds for the 
state with total magnetization 5f. = 1: 
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figure 2: The on-site magnetization < Sf > along the 
chain for different values of the biquadratic coupling, (3. 
1 was used. 



In all cases 



< V, > vbs (L) = 3 , + 2 ' 3 r , I odd. 



(4) 



1 + 3~ L 

with the general expression for any site and any even L: 

2 (l^- 1 o . y-l-L 



< Sf > vbs (L) = (-1) 



2 (ky- 1 _ o . 

-i 3 \3J 

1 - 3~ L 



(5) 



Hence, in the thermodynamic limit one obtains the fol- 
lowing result for < Sf > in the state S^ = 1: 

< sf > vbs = I {^y 1 = (-ir^e-e- 1 )^ . (6) 

This explicitly shows that the on-site magnetization is 
decreasing exponentially away from the boundary of the 
chain with a length-scale equal to the bulk-correlation 
length, £ = l/ln3 8,25 . Within the framework of a free 
boson model 26 , the same exponential decrease has been 
shown to occur at the Heisenberg point 27 , (3 — 0, in ac- 
cordance with numerical results 28,27 . Hence at the AKLT 
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point the s = 1/2 chain-boundary excitations are well- 
defined objects localized at the two ends of the chain on 
a length scale equal to the bulk correlation length. Pre- 
sumably this picture holds in most of the Haldane phase 
as we shall verify numerically in the following. 

DMRG Results: Away from the AKLT point the 
ground-state is not known and in order to study the be- 
havior of the chain-boundary excitations we use the den- 
sity matrix rcnormalization group (DMRG) method 28 . 
We refer the reader to Ref. 28 for a discussion of this 
method. The calculations were performed using the to- 
tal z-componcnt of the spin, S^, as a quantum num- 
ber and a number of states, m, between 60 and 100 
was kept. The obtained results were extrapolated to the 
m — > oo, L — > oo limit where needed. 

The bulk of our results are summarized in Fig. 2 where 
< Sf >, the on-site magnetization, is plotted as a func- 
tion of site index, i, for 8 different values of (3 between 
— 1 < /5 < 1. The calculations were in all cases per- 
formed in the Sj. = 1 state with m = 100. If the lowest- 
lying excitations are indeed determined by the s = 1/2 
chain-boundary excitations, then selecting this state cor- 
responds to polarizing the two s = 1/2 objects, and 
one would expect a rather large signal at the ends of 
the chain. This is clearly the case for most intermedi- 
ate values of (3 where well-defined excitations are visi- 
ble at the boundaries. Decreasing (3 below the disorder 
point at (3 = —1/3 clearly induces short-range incom- 
mensurate (IC) real-space correlations as the period-3 
ground-state at (3 — — 1 is approached, as can be seen 
in Fig. 2 for (3 — —2/3, —1. Increasing (3 beyond ap- 
proximately (3 ~ 0.2 moves the maximum in the on-site 
magnetization, < Sf >, away from the chain-boundaries, 
i = 1, L. As (3 approached the critical point at (3 = 1 the 
on-site magnetization resembles more and more a double- 
doughnut shape which presumably is characteristic of the 
dimerized state the system enters. 

Considering the results in Fig. 2 one may ask if 
the chain-boundary excitations remain localized at the 
boundaries as the critical points are approached. In or- 
der to obtain a partial answer to this question we have 
re-examined the results (3 = 2/3 considering a system of 
size L = 400. Our results are shown in Fig. 3. An ex- 
ponential decay is now clearly visible as detailed in the 
inset. Assuming that we can still identify the associated 
length scale with the bulk correlation length we obtain 
£ = 61(4). Hence, presumably the s = 1/2 objects re- 
main localized at the boundaries until the critical point 
(3 = 1 is reached where the correlation length diverges. 

If the s = 1/2 objects remain localized at the bound- 
aries of the chain the first-site magnetization < Sf > 
should remain finite in the thermodynamic limit. If at 
some point the s = 1/2 objects de-localize and spread 
out over the entire chain we would expect < Sf > to go 
to in the thermodynamic limit with non-trivial finite- 
size corrections. As a measure of the presence of well- 
defined chain-boundary excitations we have therefore nu- 
merically studied < Sf > as a function of (3. Two ef- 
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figure 3: The on-site magnetization in the S? = 1 state 
for a 400 site chain with (3 = |. The inset shows the same 
data for half the chain on a log-scale. The asymptotic 
exponential dependence is clearly visible, m — 100 was 
used. 

fects complicate the analysis: As (3 — > ±1 the correla- 
tion length rapidly diverges necessitating the use of larger 
and larger system sizes as well as larger values of m, and 
secondly the incommensurate correlations arising beyond 
the AKLT point. Our results for several different values 
of L, m are shown in Fig. 4. The AKLT point constitutes 
an apparent maximum in this amplitude with a value 
of < Sf >= 2/3. Where possible, we have performed 
extrapolations in both L, m to obtain results represen- 
tative of the thermodynamic limit. These are shown as 
the solid line in Fig. 4. The results are consistent with a 
non-zero value of < Sf > throughout the Haldane phase, 
vanishing as the two critical points are approached at 
(3 — ±1. The extrapolation to the thermodynamic limit 
is further complicated due to the fact that the s = 1/2 
objects strongly interact when L -C £, leading to non- 
trivial finite-size corrections to < Sf >. This is shown in 
the inset of Fig. 4 where < Sf > is plotted as a function 
of L for (3 = 2/3. A clear minimum occurs at roughly 
1-2 times the estimated correlation length. We assume 
that this minimum corresponds to a system size separat- 
ing regions with strongly and weakly interacting s = 1/2 
objects. 

Summarizing, we have presented results that show that 
s = 1/2 chain-boundary excitations occur throughout 
the Haldane phase as the biquadratic coupling is varied. 
These excitations are localized at the ends of the chain 
with a maximum in their amplitude of < Sf >= 2/3 
occurring at the AKLT point, (3 = — 1/3. It would be in- 
teresting to compare our results for < Sf >, proportional 
to the local susceptibility, with NMR results in addition 
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G- - -O After DMRG for m=1 00, L=1 00 
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figure 4: The first site magnetization, < Sf >, as a func- 
tion of (3 between — 1 < f3 < 1. The stars are data for 12 
sites (*), the circles data for L — 100, to = 100 (o), and 
the plusses extrapolations in to, L to the thermodynamic 
limit (+). The inset shows the L dependence of < Sf > 
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calculated with 



100. 



to the ESR results already available. 
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